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Abstract 

We study the entropy of Schwarzschild black holes on DGP brane. The radius of event horizon 
on DGP brane is obtained by a numerical method. It is smaller than that of Einstein gravity in 
the conventional branch, and is larger in the accelerated branch. However, the difference is very 
small. The entropy of the black hole is calculated by using the improved brick-wall method. 
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After Bekenstein has suggested that the entropy of a black hole is proportional to the 
surface area at the event horizon 1], there have been several studies for the statistical 



origin of the entropy|2|, y, IJ]. This can be understood by using a brick- wall method[5j| or 
an improved brick-wall method^, Q]. In the improved brick- wall method, a thin layer at 
the horizon is used instead of an infinite spherical box. As for a black hole on the Dvali- 
Gabadadze-Porrati(DGP) modelpl, an exact expression of metric has not been obtained. In 
this paper, we first get the radius of event horizon numerically, then, based on this data about 
horizon, calculate the entropy of the black hole by using an improved brick-wall method. 
The action in the DGP model is given by 



(1) 



S = Mi / (Px^f^gR + Mp / (Fx^-gR 



where the tilde denotes the quantities on the four-dimensional world- volume of the brane. A 
cross-over scale is defined by r c = m" 1 = Mp/(2M 3 ). We consider a spherically symmetric 
and static black hole on the brane. Then, the five- dimensional line element can be written 

as 

-e~ x dt 2 + e x dr 2 + r 2 dVt 2 + ^drdy + e a dy 2 , (2) 



di 2 



where A, 7, and a are functions only for r and y and the brane is located at y = 0. The 
equations of motion from the action ([]]) with Israel's junction conditions reduce to the 
equation for A on the brane as follows[9[]: 



U 2 - 4(1 + U)U Z - 8U{2 + U)=0, 



(3) 



where z = ln(r/r ), U = — 2 exp(— 2z)P r /?>m 2 rl } P r = dP/dr, and P(r) = r(l— exp(— A) ) | y= o 
(r being an arbitrary constant). Two solutions of Eq. ([3]), regular (or conventional) and 
accelerated branches, are given by ^ 

k 1 r = R{U) = 
k 2 r = A(U) = 



1 + 3C/ + /) 



U 2 {3 + 3U + v / 3/) 2v/5 (-5 - 3U + /). 

(-5 -3U + /)(-3 - 3*7 - V3f) 2V * 
{U + 2) 2 {l + 3U + f) 



(4) 
(5) 



where k\ and k 2 are constants of integration and / = y/1 + QU + 3U 2 . By applying boundary 
conditions by P(0) = and P(oo) = in the regular branch and P(0) = tm and P(r) ~ 
m 2 r 3 for large r in the accelerated branch, k\ and k 2 are determined as 2(r*ki) 3 = c\ ~ 0.43 
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FIG. 1: Filled boxes and unfilled boxes are the ratios of th to r^.j in the regular and the accelerated 
branches, respectively, where th is a horizon of a black hole on DGP brane and tm is a horizon of 
a Schwarzschild black hole in Einstein gravity. It shows that r# approaches tm when r c is much 
larger than tm- 

and 2(r*A;2) 3 = c 2 ~ 4.41, where r* = (rM^c) 1 ^ 3 !!^- From -P(r) = — \rr? c J dr r 2 U(r), we 



obtain 



J ] roc 

P(r)=r M 1 UR{Uf / dUR{Uf 
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(6) 
(7) 



which correspond to the regular and accelerated branches, respectively. 

Since P(rjj) = th at the horizon, we can obtain Ujj = U(rjj) from Eqs. (jlj) - 
(JU). Substituting the corresponding {7#'s into Eqs. (HJ and (jSJ), the horizons are cal- 

1 /3 

culated from r H /r M = [2(r c /r M ) 2 /ci] R(Uh) for the regular branch and r H /r M = 
[2(r c /rM) 2 /c2] 1 ^ 3 A(Uh) for the accelerated branch. The results for the regular and the ac- 
celerated branches are shown in Fig. [TJ Note that the horizon can be written as r# = Tm + 
where £ should be very small because r c j % is very large. This fact will be taken into account 
for calculating the entropy. 

In order to find the entropy of a black hole, we would rewrite the line element (J2J) on the 
brane as ds 2 = —fdt 2 + f~ 1 dr 2 + r 2 dfl 2 , where dfl denotes a solid angle. The asymptotic 
behavior of the metric is obtained as 

f(r) ~ 1 h bam 2 r 2 y—j , for r r# (8) 

/(r) « 1 - -y - 6m^r 2 , for r > r* (9) 



3 



where a ~ 0.84 and 6 = 1 and b = — 1 correspond to the conventional branch and the 



accelerated branch, respectively, 



9 



id- 



Since a = bam 2 c r 2 ^ 3_1 ^ is very small, there exists a solution r H such that f(r H ) = 0, 
where r# = + £ with |£| <C Tm- The fact that £ is very small compared to ru appears 
in Fig. CD From /(r#) = f( r M + = 0, we obtain £ « — ar^ 2v ^- Note that the radius 
of the event horizon is smaller than that of Einstein gravity in the conventional branch, 
and it is larger in the accelerated branch. Now, the Hawking temperature is given by 
T = (3~ l = k/(2tt) 7/(47rr^), where k is a surface gravity of the black hole and 7 
I - ilr H + (4 - 2V3)ar 4 H 2V * « 1 - (5 - 2V3)£/r H . 

In order to calculate the entropy of a given system in the brick-wall method, we consider 
a quantum gas of scalar particles confined within a box near the horizon of a black hole and 
introduce a cut-off parameter Q]. The free scalar field is assumed to satisfy the Klein-Gordon 
equation, (□ + m 2 )\I> = 0, with boundary conditions ^{rn + h) = ^(L) = 0, where Th is the 
horizon, m is the mass of a scalar field, h is an infinitesimal cut-off parameter, and r# + h 
and L represent the inner and the outer walls of a "spherical" box, respectively. Suppose 
that this system is in thermal equilibrium at a temperature T with an external reservoir. 
Using \l/ = exp(— iut)tp(r, 6, <p), the Klein-Gordon equation is reduced to 



d 2 i\) {2 1 df\ dijj 1 



dr 2 \r f dr J dr ^ f 



2 1 f & nd 2n d 2 

— - m 2 + — — - + cot 6— + esc 2 9 



V> = 0. (10) 



/ r 2 \d6 2 06 dcp 2 i 

By the WKB approximation, using ip ~ exp[iS(r,9,<p)], we obtain p 2 = f~ l [oj 2 /f — m 2 — 
pl/r 2 — p^p/ir 2 sin 2 #)], where p r = dS/dr, pg = dS/d9, and p^ = dS/dip. Then, p 2 = 
g rr p 2 + g 06 Pe+ 9 w vt> = oj 2 / f — m 2 . The number of quantum states with energy not exceeding 
uj can be written as 

If 2 f r 2 ( J 2 \ 5 

n(cu) = J drdOdipdprdpgdp^ = — J dr — I — - m 2 \ . (11) 

The free energy is given by 

r , n(u) 2 f , r 2 f°° , (x 2 -x 2 ) 3/2 
F = - du a y ' = — / dr— / dx± ^ — , 12 

where x = [3uj and xq = (3ra\ff '. 

The degrees of freedom are dominant near horizon. Thus, we consider a thin layer with 
L = th + h + 8, where S is small. Since xq goes to zero near horizon, we obtain 

2 r , r 2 f°° , x 3 2tt 3 f rH+h+s , r 2 



3tt/? 4 y - ; 2 7 & - 1 45/? 4 y rH+/l p 



The metric near the horizon can be written as /(r) ~ 7(1 — Th/t). Therefore, the free 
energy can be calculated as 

27r 3 ri 5 . lS 

F ps (14) 

45/? 4 7 2 h(h + ti)' y } 

and the entropy becomes 

S =f? Q JLn *!*JL I (15) 

1 8(3 45/3 2 7 2 /i 2 (/i 2 + 5 2 )' 1 ' 

where h = J^ H+h ^/g^dr, 5 = J^ H ^ +S \Jlh~rd' r - Note that h and 5 are the cut-off parameter 
and the thickness of the layer, respectively. Thus, we can write the entropy as 

1 1 S 2 

s ~i A '9o^h 2 (h 2 + py (16) 

where A = Anrfj. When the cut-off parameter h is chosen as S 2 /jh 2 (h 2 + 5 2 )] = 90tt/G, Eq. 



T6|) agrees with the Bekenstein-Hawking entropy 5bh = 
One may consider the cut-off introduced in the brick-wall method a bit ad hoc. In Ref. 



11 



121 ]. this point was criticized and the Pauli-Villars regularization scheme was used to replace 



the cut-off. In Refs. 13| . the entropy was evaluated based on the notion of entanglement 
entropy. Below, we evaluate the entropy in the thin layer using the Pauli-Villars scheme and 
compare it with our result (Tl6i) . Introducing the five regulator fields (i — 1, ... ,5) besides 
the original field (i = 0) [11] , the total free energy becomes 

F«-pL± Al r _^ f d y ~^f\ (17) 

3l i Jo "'"-'i S 2 (l-S) 4 

where A = A 3 = A 4 = +1 for the commuting fields and Ai = A 2 = A 5 = — 1 for the 
anticommuting fields and the integration is taken for values where the square root is real. 
Setting hi — and focusing only on the divergent contributions at the horizon, we obtain 
the entropy 

S ~ ^ ( B I ^ a] - ABl I Al (18) 
b ~ 3(3 [ B+ 15(3^)^ 4 12vr + 90' 



with the same A and B coefficients given in Ref. [11]. As it was discussed in Ref.jll|, the 
second term including A corresponds to the 1-loop correction to the quadratic-curvatre 
terms of the gravitational acton, and is irrelevant to our case. The first term corresponds 
to the 1-loop renormalization of the Bekenstein-Hawking entropy S = A/AGr, where the 
renormalized Newton's constant Gr is given by 1/Gr = 1/Gb + B^/12tt with the bare 
Newton's constant Gb- Note that our result slightly differs from the result of Ref. 
factor 7. 
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